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The behavior of a neutral polyampholyte (PA) chain with A'^ monomers, in an ionic solution, is 
analyzed in the framework of the full Debye-Hitckel-Bjerrum-Flory {DH BjF) theory. A PA chain, 
that in addition to the neutral monomers, also contains an equal number of positively and negatively 
charged monomers, is dissolved in an ionic solution. For high concentrations of salt and at high 
temperatures, the PA exists in an extended state. As the temperature is decreased, the electrostatic 
energy becomes more relevant and at a T = Te the system collapses into a dilute globular state, or 
microelectrolyte. This state contains a concentration of salt higher than the surrounding medium. 
As the temperature is decreased even further, association between the monomers of the polymer 
and the ions of the salt becomes relevant and there is a crossover from this globular state to a 
low temperature extended state. For low densities of salt, the system is collapsed for almost all 
temperatures and exhibits a first-order phase transition to an extended state at an unphysical low 
temperature. 

PACS numbers: 36.20.Ey, 64.60. Cn, 82.35.+t, 87.15.By 
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I. INTRODUCTION 

In recent years, there has been a lot of interest in 
the study of the physical properties of macroniolecules. 
When dissolved in a good solvent, a neutral macro- 
molecule, that is made of A'^ units (monomers), is char- 
acterized by a radius of gyration that scales as i? ^ N^a, 
with ly = 3/5. In a bad solvent, as the temperature is 
lowered, a polymer undergoes a phase transition from an 
extended state, characterized hy — 3/5, to a compact 
globular state, with i/ = 1/3. At the transition point, 
T = Tg, the system assumes a coil configuration, with 
V — 1/2. The transition resembles the tricritical behav- 
ior present in some magnetic materials 01 and can be 
described by the Flory-De Gennes theory^-^ . 

In nature, however, charged polymers are prevalent. 
They can be divided into two classes, polyelectrolytes 
and polyampholytes {PAs) . The most common example 
of a polyelectrolyte is a DNA dissolved in water. Upon 
dilution, the phosphate groups of the DNA molecule lose 
their atoms and become ionized. All the ionized 
groups of the polyion carry the charge of the same sign. 
The repulsion between the like charges helps to stretch 
the molecule and contributes to the overall persistence 
length of a polyion |||,^ . 

Proteins, on the other hand, are an example of polyam- 
pholytes. A PA molecule can have monomers of either 
positive or negative charge. Thus, two extremes are pos- 
sible. In one extreme, we can have a neutral PA, in which 
the number of positive monomers is exactly equal to the 
number of negative monomers. In the other extreme, all 
the monomers can be of the same sign, in which case a 
PA will be equal to a polyelectrolyte. In this paper we 
shall confine our attention to the case of a neutral PA. 

In contrast to a normal polymers, the understanding 



of PAs is rather poor. This is due to the difficulty of 
properly accounting for the long range Coulomb inter- 
action between the monomers of a PA. The long range 
electrostatic interaction is qualitatively different from the 
short-range excluded volume repulsion, that drives the 
coil-globule transition in an usual uncharged polymer. 
Consequently, our understanding of the full conforma- 
tional structure of the PA is rather incomplete. 

One of the first attempts to study a neutral PA was 
done by Edwards et al. 0. They argued, on the ba- 
sis of the Debye-Hitckel {DH) limiting law, that when 
the polymer chain is in a nonionic solvent, the competi- 
tion between the net attractive electrostatic energy of the 
PA and the reduction in its entropy, due to the confine- 
ment in a sphere smaller than its natural radius, leads 
to a transition from an extended coil to a "dilute glob- 
ular state" (microelectrolyte). The microelectrolyte is 
characterized by a density smaller than that of the col- 
lapsed state of a chain made of neutral monomers. The 
microelectrolyte picture of a PA was further extended 
and elaborated in the work of Higgs and Joanny |^ , who 
have concluded that a sufficiently long PA will always be 
collapsed. 

In a recents papers, we have investigated the pos- 
sibility of a different type of coil-globule transition in a 
neutral non- alternating PA. In the absence of salt, the 
behavior of a neutral PA, with N charged monomers, 
was analyzed in the framework of the full DH theory 
[|lH , augmented by Bjerrum's idea of ion association |l^] 
and Flory's affine network theory of rubber elastic- 
ity. In this approach, the allowance was made for the 
fact that a strong electrostatic attraction between the 
unlike monomers makes it favorable for them to associate 
forming intermolecular bridges. The PA molecule then 
resembles an affine network with the crosslinks of func- 
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tionality 4. It was found that at high temperatures the 
crosshnks are quite rare, and the picture of a PA as ex- 
isting in a dilute globular state remains unchanged, with 
the radius of gyration scaling as R ^ aN^ and v ~ 1/3 
(microelectrolyte) , as predicted by Edwards et al. [0,1). 
As the temperature is decreased, the ionic association 
becomes energetically favorable. At very low temper- 
atures, the fraction of unassociated charged monomers 
approaches zero and the PA molecule resembles a mi- 
crogel, with v = 2/b (compare with the Flory exponent 
Vp = 3/5). The two states are separated by a first order 
phase transition. It is interesting to note that unlike the 
case of an usual neutral polymer, the extended state of 
PA [v = 2/5) occurs at low temperatures. 

In order to analyze the effects of addition of salt to 
the polymeric solution, Higgs and Joanny |^ extended 
Edwards et al. j^] method, including the electrostatic in- 
teraction between the charged monomers of the chain and 
the free ions of the salt. They considered the case where 
the density of charges on the PA is of the same order 
of magnitude as the density of added salt. They have 
concluded that the only effect of added salt was to renor- 
malize the second virial coefficient, adding an attractive 
term to it. At high temperatures and in a good solvent, 
the excluded volume term dominates the attractive DH 
contribution and the chain is extended. As the tempera- 
ture is decreased, the electrostatic contribution balances 
the excluded volume contribution and the polymer col- 
lapses into a dilute globular state. 

Higgs and Joanny's theory (HJ), however, is based on 
the assumption that the density of the polymer is close 
to the density of salt, what is actually a good approxi- 
mation only when the densities are quite high. Further- 
more, since HJ work in the framework of the linearized 
DH theory, their conclusions can only be trusted at high 
temperatures and cannot be extended to the low temper- 
atures, where they find an unphysical divergence of the 
concentration of ions near the PA. 

In order to further explore these issues, we analyze 
this problem in the framework of the full Debye-Hiickel 
Bjcrrum-Flory (DHBjF). The possible association be- 
tween the counterions and the monomers of the chain 
is explicitly taken into account. In a good solvent and 
for density of salt po > p (p is the density of polymer 
in extended state at high temperatures), corresponding 
to the region where HJ found the transition, the system 
has the following behavior. For high temperatures, the 
repulsion due to the excluded volume dominates over the 
DH electrostatic attraction and the polymer is extended. 
The counterion association is insignificant in this region. 
As the temperature is decreased, the electrostatic energy 
prevails and at T = Tg the system collapses into the di- 
lute globular state, in agreement with HJ predictions. 
We obtain explicitly the values of the transition temper- 
ature for different densities of salt and for different vol- 
ume fractions of charged monomers. As the temperature 
is lowered even further, the counterion-monomer associ- 
ation becomes energetically favorable. The number of 



ions of salt attached to the chain increases. As T ^ 0, 
the fraction of unpaired monomers approaches zero and 
the PA molecule becomes extended like an usual neutral 
polymer. Within our approach we find that the instabil- 
ity observed by HJ is hidden inside the metastable region 
of the salt solution. Furthermore, beyond the HJ results, 
we are able to investigate the properties of a polymer- 
salt solution when the density of salt is quite small. We 
find that for po < p, the fluctuation induced attraction 
between the monomers is strong and the polymer exists 
in a dilute collapsed state for almost all temperatures. 
At very low temperatures and low densities, the system 
undergoes a first-order transition into an extended state. 

The paper is organized as follows. In Sec.|| our model 
for PA-|-salt is presented and the Helmoltz free energy 
of the composed system is constructed. A full thermo- 



dynamic analysis is formulated in Sec, III, in particular, 
the equations that determines the equilibrium configura- 
tions of the system. In Sec. IV a detailed discussion of 
our results is presented. 



II. DESCRIPTION OF THE MODEL 



A. Introduction 



We consider a polyampholyte {PA) chain consist- 
ing of N constitutional units (hereafter referred to 
as monomers) of three different types: iVjv neutral 
monomers, {N — Nn)/2 positive monomers and {N — 
Nn)/2 negative monomers, randomly distributed along 
the chain, each spherical in shape with a diameter a. 
The polymer chain occupies an effective sphere of vol- 
ume Vi = AttR^ /3, where R, the radius of gyration, mea- 
sures the end-to-end extension of the polymer. The re- 
duced density of monomers, p* = pa^, inside the globule 
is p* = 3A^aV(47ri?3). 

The solvent that occupies the total volume V of the 
system is represented by an uniform medium of a di- 
electric constant e in which an amount of simple 1:1 
salt has been added. The positive and negative ions 
of salt dissolved in the solvent have a fixed densities 
Po ^ Po ^ Po ^ "^Po —2/3(7. The basic approximations 
of our model are: 

1. In the limit of large dilution wc can neglect the 
interactions between different PAs and concentrate our 
attention on one chain. 

2. The most relevant interaction between the 
monomers of the chain and the free ions of salt is of 
electrostatic nature. For simplicity, the excluded vol- 
ume interactions between the PA and the salt will be 
neglected. 

3. The short range interactions between the polymer 
and nonionic solvent are included in the effective value 
of the second virial coefficient. 
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FIG. 1. Schematic drawings of a PA chain with neu- 
tral monomers O) (-^ ~ Nn)/2 positive monomers ©, and 
(A*' — Nn)/2 monomers 0, immersed in an ionic solution. 
The salt ions can associate with the positive and negative 
monomers of the chain, forming dipolar pairs. The PA is 
confined to a spherical region of radius R. 

The strong electrostatic interactions make it favorable 
for the ions of salt to associate with the PA. After 
the thermal equilibrium is achieved, the chain is com- 
posed of three different structures: neutral monomers, 
free charged monomers and monomer that have an asso- 
ciated ion of salt attached to them (see Fig.|l]). We shall 
refer to the latter as a "dipolar pair" . The total number 
of monomers can then be written as 

N = NN + N++N-+Nd, (1) 

where iV+ and N~ are the number of unassociated posi- 
tive and negative monomers and Nd is the total number 
of dipoles. Defining / as the number fraction of charged 
monomers in the PA, and x as the fraction of charged 
monomers that form dipoles, we can write 

ATjv = (1 - f)N , neutral, 

N+ = |(1 - x)N = N~ , charged, 
Nd = fxN , dipoles. (2) 

It is then possible to express the densities of neutrals, 
charges and dipolar pairs, in terms of the total density 
of the chain, p = N/Vi, by the expressions 

PAT = (1 - f )p , 
f 

Pd = fxp ^ Pd + Pd ' (3) 

where p'^ and are the densities of dipoles containing 
positive and negative monomers, respectively. In prin- 
ciple, ions of salt are different from the monomers of 
the PA and the dipoles formed of positive and negative 



monomers have to be distinguished. Since, however, the 
PA is neutral, in equilibrium Pd = Pd ■ 

The salt particles are free to move along the whole 
volume of the composed system. In the absence of the 
PA, the positive and negative charges are uniformly dis- 
tributed with a constant density po. In the presence of 
PA, due to electrostatic fluctuations, the density of salt 
near the the polymer is not equal to po- Therefore, in or- 
der to treat the problem in a simple way, we shall divide 
the total system of volume V into two subregions: "in- 
side" and "outside" the volume Vi. The density of salt 
particles outside is uniform and equal to po = Nq/Vq, 
where V = Vi + Vo. The density of the salt inside the vol- 
ume Vi is ps = pt + Pj ■ Since there is an equal number of 
positive and negative free ions, pf = pj = ps/2 = Ns/Vi. 



B. Helmoltz free energy 



The Helmoltz free energy for the whole system is the 
sum of two terms, F — Fi + Fq. The "inside" free energy 
Fi , associated with the subregion where the polymer and 
the salt coexist, and the "outside" contribution Fq, for 
the regions where only salt is present. The free energy in 
the outside region is composed of: 
. an energy due to short range repulsion; 
. an electrostatic energy due to the interaction between 
ions; 

. and the kinetic-entropic contribution. 
For the inside region, the Helmoltz free energy Fi, con- 
tains: 

. an elastic contribution due to the elongation and con- 
traction of the chain; 

. an energy due to the monomer-monomer interaction; 
. an energy due to the ion-ion short range repulsion; 
. an electrostatic energy of monomer-monomer, salt-salt 
and monomer-salt interactions. 

. an entropic energy due to the mixing of the different 
types of particles present inside the region Vi . 

a. Elastic free energy 

According to the Flory-De Genncs [p|, p^|jl^ ] theory, the 
elastic free energy of the chain is 

/5Fe/ = ^(a2-l)-31na, (4) 

where a = R/Rq is the expansion factor of the chain, 
measured relative to the non-strained Gaussian state 
with radius Rq ~ a\/N . 

b. Free energy of excluded volume 

Following the usual Flory-De Gennes theory, the 
short range interactions between the particles of the 



3 



chain, at low densities, are approximated by a virial ex- 
pansion 



(5) 



where Wi and W2 are second and the third virial coeffi- 
cients. 

The excluded volume interaction for the salt is approx- 
imated by the free- volume form 



PFlic = ~V^ps lii{l-Bps 



(6) 



where B is chosen to yield appropriate maximal pack- 
ing densities or to match high-T second virial coefficient. 
Since we are assuming that the relevant interaction be- 
tween the monomers and the salt particles is electrostatic, 
we do not consider the excluded volume repulsion be- 
tween the free ions inside Vi and the monomers. Outside 
Vi, the excluded volume contribution to the Helmoltz free 
energy is 



HC 



(7) 



and, for small densities po, it can be neglected. 

c. Electrostatic term 

The electrostatic interaction between the charges is 
calculated in the framework of the Debye-H-itckel {DH) 
theory. The standard argument leads to pT| , p^| -p8[ 



T*[Ka) 



ln(Ka + 1) 



{nay 



(8) 



Inside Vi, the inverse Debye screening length k = 
is determined by all free unassociated charged species, 
i.e., both monomers and salt ions, so that (na)^ = 
4np*jT*, where pi = a^Ni/V = p* + /(I - x)p* 
and T* — ksTea/q^. Outside the volume Vi, the De- 
bye length is determined by the concentration of salt, 
{naf = {noaf = Anpl/T*. 

The DHBj theory has proven to be successful in giv- 
ing the correct qualitative picture of a simple electrolyte 
pq-ESl. Thus, while inclusion of higher order contri- 
butions, such as dipolc-ion and dipole-dipole interac- 
tions, improves the quantitative agreement of the th eory 
with experiments and Monte Carlo simulations pO[ , 
most of the important physics is already captured at the 
DHBj level |L7| . As a leading order approximation we 
shall, therefore, treat the dipolar clusters as ideal nonin- 
teracting specie. 

d. Free energy of mixing 

Inside the volume Vi , the system is composed of a mix- 
ture of neutral, positive and negative monomers, as well 
as, positive and negative ions of salt, and dipoles made of 
monomers and associated counterion. If the bonds con- 
necting the monomers were cut, we would have a "gas" 



made of these species, 
mixture is 



The ideal gas free energy of the 



p,K 



(9) 



where A.j{T) is the mean thermal wavelength, Aj{T) = 
h I ^Ji'K'fhjkBT , fhj is the geometric mean mass for a clus- 
ter of kj units and (T) is the internal partition function. 
The index j specifies all the possible species; positive, 
negative, and neutral monomers, dipoles, and positive 
and negative free ions. In our case, fc^r = fc- = fc+ = 
ks = 1 and ^2 = 2. To simplify the notation, we shall 
take the masses of both monomers and ions to be equal 
to m; therefore, Aj(T) = A. Since all the monomers 
are hard spheres, the molecular partition is a constant, 
= £,1 = £,s = £, ^ ^- For dipoles, however, the internal 
partition function is |I2] 



^2(r*) = 47r/ e'^^'^'r^dr = K{T*) 

J a 



(10) 



where K{T) defines an association constant for the for- 
mation of bound dipolar pairs. For the cutoff c, we use 
the Bjerrum's choice, a/2T*, which corresponds to the 
inflection point of the integral as a function of c. The 
Bjerrum equilibrium constant is 

Ei \ — ] - Ei{2) + 



K{T* 



3(T*)3 _ 



l/T* 



1 



1 



{T*y 



(11) 



As was shown by Ebeling the low temperature 
asymptotic value of this constant is in exact agreement 
with the cluster expansion . 

While the ions of the salt are free to move indepen- 
dently of each other, the monomers of PA are constrained 
by the bonds that maintain the integrity of the chain. 
The equation (|^), therefore, overestimate the entropic 
free energy by the amount equivalent to the the free en- 
ergy of the ideal gas composed of N particles. The free 
energy of mixing is therefore. 



MIX — 



E 



iV.ln 



p,A, 



N, 



= VpN In A>Ar + Vp+ In A>+ + Vp 

Pd 



InA^ 



-l^,p+lnA" 



Pd 



Vp^ In A" 



K{T) 



InA^p; 



KiT) 

+Vpt \nA^pt + Vp; 
-VplnA^p - Vpt ~ Vp- . (12) 

Outside, in Vq, the entropic contribution to the free 
energy is given by an ideal gas free energy 

= Vop+ In A3p+ + Vopo In A^po - Vop+ - K)Po ■ 

(13) 
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III. THERMODYNAMICS OF THE SYSTEM 



Mo = Ms 



(18) 



The equilibrium configuration of the P^+salt system 
is determined by the minimization of the Gibbs free en- 
ergy. Since the total system is composed of two sub- 
regions, the total Gibbs free energy can be written as 
G = Gi-t-Go, where d is the Gibbs free energy of the sub- 
system formed by the "inside" region and Go is the Gibbs 
free energy for "outside" region, Vq. The minimization 
of the total free energy imphes that 5G = SGi + 6Go — 
and, consequently. 



= 



dG, 



dG 



dV,. 



dGj 
'dVo 



dG, 

'dm 



-dNj 



dN: 



-dVo- 


dG, 


dN+- 


dG, 


dN+ 


dNp 


dG, 


dN^^ 


^ dG, 


dN+ 


dN^ 


dN+ 


dGo 


dN+ H 


dGo 


dN, . 


dN+ 


dNo 



dN- 



(14) 



The total volume of the system is fixed, V = Vi + Vq, 
which leads to the condition SVi — —SVq and, conse- 
quently, to the equality of the osmotic pressures P = 
—dG/dV inside and outside the globule, 



P, = P„ + Ps = Po 



(15) 



In the above expression, our first condition for equilib- 
rium, Pp,Ps and Po are the osmotic pressures of the poly- 
mer, salt inside the volume Vi, and salt in Vq, respec- 
tively. 

Now, since the total number of ions of salt, N~ + 
^d+^d+ ^0 + ^0 ' and the total number of monomers 

+N++Nj;+N+ (inside Vi), are 
-SNd (where Nd = N^^ + N~ 



that contain charges, 
constants, we have i5iVr, = 



and Np = N' 



N+) and SNo = SNd - SNs (where 
Ns = N;r + N+ and No = N, + iVo+). Furthermore, 
since we are considering a neutral PA and a neutral salt 
solution, the chemical potentials, fif = dG/dNf and 
= dG/dNp , for negative and positive species are 
equal, /i^ — fi'f — and /i~ = /i+ = /ip. Therefore, 
Eq.([l4|) becomes 

= -tipdNd + tXddNd + ^lsdN, - no{dNs + dNd) = , 

(16) 

where Mp and Md are the chemical potentials of the 
charged monomers and the dipoles along the chain and /is 
and Mo are the chemical potentials of free salt inside and 
outside the volume Vi, respectively. Since the number of 
dipole pairs and free ions are independent variables, the 
Eq.(|l6|) leads to two additional equilibrium conditions: 
the chemical association. 



Mp + Ms = Md 



(17) 



and the condition for free exchange of salt between the 
"inside" and the "outside" regions. 



Instead of working with a full system composed of both 
the "inside" and the "outside" regions, we can confine our 
attention strictly to the "inside" region, at the expanse 
of introducing some Lagrange multipliers. This leads to 
an effective Gibbs free energy 



f3g = f3F + (3PoV - f3fioNs - (3fioNd 



(19) 



where F is the Helmoltz free energy of the "inside" re- 
gion, 

(3F = f3Fei + iSF^c + PFfjc + PFdh + PFmix , (20) 

and the pressure Po and the chemical potential mo play 
the role of Lagrange multipliers. All the thermodynamic 
information about the "inside" region is now contained 
in g. The equilibrium state of the PA + salt subsystem 
is determined by the minimization of G, 6Q = 



d(3F 
dV 



+ m ] dV 
(d(3F 



/ dpF 
[dNs 



/3mo dNs 



\dNd 



/?Mo dNd 



dN-/''^ 







(21) 



Making once again use of the conservation of the number 
of charged monomers and of the number of particles of 
salt, it is not hard to show that this reduce to three 
equilibrium equations (|l|), (0) and (^, respectively 

To find the quilibrium state of the P^-|-salt system, 
these equations must be solved simultaneously. Explic- 
itly, the conditions for equality of pressures and for free 
exchange of salt between the "inside" and the "outside" 
regions are respectively 



3(a2 



-I- 



1 



87ra3 
1) , 



21n(z + 1) - 



ziz + 2) 



1 



Bpl 



47ra3a3iV3/2 
1 



2 1 - Bpo 



87ra3 



21n(zo + 1) 



zo{zo + 2) 



1 



(22) 



where z = ku — ^ At:p\/T* and zq — kqa = A'Kp'^/T* 
and, 



(1 - Bps)po 
(1 - Bpo )ps 



exp 



1 



Zo 



z+l 
Bpo 



Zo 



Bps 



I -Bpo 



Bps 
(23) 



The condition of chemical equilibrium between the 
dipoles and the monopoles reduces to. 



2x 



KiT*) 
(1 - x)ps 1 - Bp, 



■ exp 



1 



Bps 



T* z + l l-Bps 



(24) 



The equilibrium configuration of the PA, for a fixed T 
and Po, is given by the solutions of the coupled Eqs.(p2[), 
(H) and M). 
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IV. RESULTS AND DISCUSSION 



For a fixed temperature and the external density po, 
these equations can now be solved to yield the values of a, 
X and ps- However, as the temperature is lowered, below 
T* (see Appendix^ for details), the salt solution phase 
separates into a high (liquid) and low (vapor) density 
phases. A PA then finds itself either in the high or low 
density phase of the salt solution (sec Fig.^). We shall 
assume that the PA will be in the high density phase if 
the original density pp > Pc 7 where is the critical point 
density for the pure salt solution ||l6|,|l^, and in the low 
density phase pg < p*. We chose a path such that, as the 
temperature is lowered, the density of the salt remains 
constant and equal to po up to T* (the temperature at 
which Po — pcy{T*), where pa- is the coexistence curve 
boundary). Below T;, po = Pa{T*). 

In order to obtain the phase-diagram of the system, 
we assumed Wi — Aira'^/S, the value corresponding to 
the second virial coefficient for the gas of hard spheres, 
A'' = 100 and / = 1. We chose two different paths: one 
on the vapor side of the coexistence curve, represented 
by the path V in Fig.^ and another on the liquid side of 
the coexistence curve, represented by the path L. 



ated with a = a, where a is the solution of Eq.(|27|), with 
Wi> (see FigJI(b)). Solving Eqs.(|2|), 0) and (|2|), 
the following phase-diagram arises (see Fig]^and Fig]4). 
At high temperatures, the excluded volume repulsion is 
the dominant interaction and the polymer is extended 
(scaling as = 3/5). As the temperature is decreased, 
the fluctuation induced electrostatic attraction increases 
and, at a temperature Tg the polymer collapses into a 
dilute globular state, the microelectrolyte (see Appendix 
^ for details). As the temperature is decreased even fur- 
ther, association between the monomer and the counteri- 
ons begins to takes place. The values of x and a smoothly 
increases and the polymer expands. Meanwhile, the den- 
sity of the salt inside the globule approaches from above 
the value of the density outside the volume Vi, p\{T). 
In the Appendix ^ this collapse transition is analyzed in 
detail and the transition temperature Tg* is calculated in 
the framework of the full DH theory for different values 
of salt density and fraction of charged monomers. The 
unphysical divergence of the density of salt in the vicinity 
of a PA, encountered by HJ is no longer present, being 
hidden inside the coexistence region, as shown in Fig.^. 

^ V EXTENDED COLLAPSED 

(a) I \ ^ 

t' 




FIG. 2. The coexistence curve for the simple 1:1 salt, 
for the restricted primitive model (RPM) according to 
Debye-Hiickel (DH) theory, with hard-core repulsive contri- 
butions. The dotted box is represented in the inset. The crit- 
ical point (open diamond) is T* « 0.0624 and p* ~ 0.00462, 
which implies two paths for the density of salt outside Vi: 
{V) on the vapor side of coexistence curve (open circles) and 
(L) on the liquid side of coexistence curve (open squares). 
The dashed line represents the HJ line of instability, and the 
dot-dashed line is the minimum density p of the polymer. The 
solid line represents the 9 temperature Tg. 

First, the high temperature density of the salt outside 
the volume Vi was fixed at po = 0.0139 > p (correspond- 
ing to T* — 0.06), which implies a quench on the liquid 
side of the coexistence curve as shown in Fig.^. Here, p is 
the density of the polymer in the extended state, associ- 



/•t V DILUTED COLLAPSED EXTENDED 

(b) I , . 

i;* T* 

FIG. 3. Schematic phase-diagram for: (a) po on the vapor 
side and (b) po on the liquid side of the coexistence curve. 
The temperature is very low. 

Next, the high temperature density of salt outside Vi 
was fixed at pj = 0.003 < p* (at T; = 0.0618), which im- 
plies a path on the vapor side of the coexistence curve. 
The solution of Eqs.(p^), ( p3| ) and ( p^ ) leads to the phase- 
diagram for this region, shown in Fig.^ and Fig.^. Even 
for high temperatures, the electrostatic energy is quite 
large and, consequently, the polymer is found in the di- 
lute globular state for any physical temperature. There- 
fore, no high temperature collapse transition is found, 
when the density of salt is pg < p. As the temperature is 
lowered, however, the ionic association becomes energet- 
ically favorable and the free ions of the salt to condense 
onto the charged monomers of the PA. Indeed, in this 
region x, ps and a are not a single valued function. Be- 
sides the dilute globular state (designed with d in Fig||), 
there are two other branches: one metastable (designed 
with m in Fig.^ and another stable branch (assigned by 
I in Fig.||). In the presence of a good solvent, the poly- 
mer in this last stable branch is actually extended. At 
T* TuT* (see FigJI), there is a first-order phase transition 
between the dilute globular state and the low tempera- 
ture extended state. This transition actually happens at 
such a low temperature (see Appendix ^ for details), that 
it is even difficult to give a good estimate numerically. 
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FIG. 4. (a) Fraction of dipoles x, (b) expansion factor 
a and (c) density of salt inside of spherical region, ps, for 
a quench on the liquid side of the coexistence curve, is 
shown against the temperature, for = 100, f — 1 and 

For a quench on vapor side, as the temperature is de- 
creased, the density of the salt inside Vi, Ps, exhibits 
an interesting behavior. Even if the density outside Vi 
is fixed at the vapor phase, the density inside Vi is at 



p\{T*) (see Fig.p(c)). The salt actually separates in two 
phases, a liquid phase inside the globule and a vapor 
phase outside. Therefore, even for a very low density of 
the salt outside the globule, po = the density of 

salt inside is quite large {ps = p^^{T*)). 

Finally, the density of the salt was fixed at pi < Pq = 
0.00698 < p (at T* = 0.0618), which imphes a quench 
on the liquid side of the coexistence curve, but in the 
region where the collapse transition is absent. In this 
case, like in the previous situation, the polymer is in the 
dilute globular state for almost all temperatures. Differ- 
ently from what has happened on the vapor side of the 
coexistence curve, at low temperatures there is a smooth 
crossover from the globular state to the extended state. 

As we pointed out in the last two cases, when the den- 
sity of the salt is small, po < p, the DH contribution is 
governed by the density p of polymer. For large tempera- 
tures, there are almost no dipoles, a; — s- and the density 
of salt inside the volume Vi approaches the value outside 
Po- In this region z becomes small and the DH limiting 
law can be used in Eq.(p2|), which becomes 



3(a2 



1) 



47ra3a3iV3/2 



+ Wep'/^ = , 



(25) 



where We = fa^''^^^f^/Q{T*f/'^. In the iV ^ oo limit, 
the solution of Eq. (Eq) implies a finite value for the den- 
sity p and, consequently, the size of the polymer scales 
as 



R = a^N^I^ , 



(26) 



where ao = iT* [W^ /Att) 'a/ J{1 - x), with = 
Wi/a^ as a dimensionless parameter, corresponding to 
the quality of the solvent. 

Bellow the conformation of the PA is governed by 
branch I (see Appendix |c|) . Since along this branch all 
the monomers are neutralized by the counterion, the at- 
tractive electrostatic interaction is absent and in Eq. ( |2^ ) 
the hard-core repulsive term (repulsive) has to be bal- 
anced by the elastic term (attractive). In this case, we 
find that 



3(a2 _ 1) 



W2p-^ = 



47ra3a3iV3/2 
which is the same as for a neutral polymer: 



(27) 



(i) If Wi 3> (good solvent), the elastic term — 3(a^ — 
l)/47ra3a3Af3/2 balances the hard-core contribu- 
tion Wip^/2 and gives a radius of Gyration R — 
aiN^/^, where the amplitude oi = {3W^ /Stt)'^^'^ a. 
So, the PA is in an extended state; 

(ii) if Wi ~ (theta solvent), the second virial term is 
zero, and the elastic term balances the third virial 
W^p^ ■ In this case the radius scales as i? = a2iV^/^, 

where as = {%W;^ j {A-nff^^ a, and 1^2* = W^2/a^ 
In this case, the PA behaves like a coil; 
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FIG. 5. (a) Fraction of dipoles x, (b) expansion factor 
a and (c) density of salt inside of spherical region, pa, for 
a quench on the vapor side of the coexistence curve, is 
shown against the temperature, for = 100, f ~ 1 and 
Wi = 4wa^/3. The stable branch d represents a dilute col- 
lapsed state, I represents a stable extended state and m is a 
motastablc state. The dashed line in (c) represents the coex- 
istence boundary. Notice how the density of salt inside the 
globule asymptotes to the liquid side of the phase boundary. 



(iii) if Wi ^ (poor solvent), the second virial coeffi- 
cient \Wi\p'^/2 balances the third virial W2p^ and 
the radius of Gyration scales as i? = asN^^^, where 
03 = {3Wi/2Tr\W^\f^^a. Thus, the PA is in a 
dense collapsed state. 

Unfortunately this interesting behavior will be impossible 
to observe experimentally since is too low. 



V. CONCLUSION 



In this work, wc have analyzed the behavior of a neu- 
tral polyampholyte in the presence of salt, using the full 
Debye-Hitckel-Bjerrum-Flory (DHBjF) theory. Allow- 
ing for the association between the charged monomers 
and the free ions of salt, we found that when the density 
of salt is large po > p, at high temperatures the PA is 
extended. As the temperature is decreased, there is a col- 
lapsed transition to a dilute collapsed state at T = Tg. As 
the tempcratiuc is decreased even further, the free ions of 
the salt associate with the charged monomers, decreasing 
the electrostatic attraction. In fact, at zero temperature, 
all of the monomers are neutralized and, since there is 
no free charges left, the chain behaves as a usual neutral 
polymer. Therefore, in the presence of a good solvent, 
there is a crossover from a dilute globular state to an ex- 
tended state at low temperatures. If the density of salt is 
small, po < p, the polymer is in the dilute collapsed state 
for almost all temperatures. A first-order phase transi- 
tion to a low temperature extended state occurs for an 
unphysical low temperature. 
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APPENDIX A: COEXISTENCE CURVE FOR 
THE SALT 

The salt solution exhibits two possible phases, liquid 
(high density) and vapor (low density). The coexistence 
curve boundary Pa{T*) between these phases is given by 
the solutions of the coupled equations 

P{pl) = P{p',), (Al) 

and 

M=M, (A2) 
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where P = -dpP/dV and fi = d(3F/dN. The Hel- 
moltz free energy /3F = PFhc + fiFoH + PPid is given 
by the Eqs.(0), (|) and Exphcitly, 



Bpl 



1 - Bp„ 87ra3 



21n(z + 1) 



z + 1 
(A3) 



and 



P{P<y) 



1 



2r* z + 1 

I -Bp,' 



InA^^ 



ln(l - Bp,) 



(A4) 



where z = A/47r/?jyT*. The resuhing coexistence curve, 
for B = Aa^/Sy/S (bcc packing), is shown in Fig|| and, 
for low temperatures, the vapor and hquid branches are 
given analytically by 



Pa ~ exp 



and 



1 



1 



1 



Pa ~ Pm 



2T* 2 V ttpmT* V 2 SnpM 



PM 8ttpm 



1 - W^PmT* - 2T* In^ 
^{2lnAnpM-2-8npM)T* 



InT* 



(A5) 



(A6) 



where pM = l/B is the maximum packing density. As 
the temperature T* T* , p'', - p^ ^ and at the 
critical point, p"^ = p^^ — p*. The critical temperature 
and density are T* « 0.0624 and p* « 0.00462. In this 
case, the reduced temperature is only 10% higher than 
the MC estimate (T* = 0.057 ± I5, p* = 0.030 T 8) [|o|, 
but the critical density is too small. The low value of 
critical density is due to the fact that DH theory omits 
important nonlinear configurations, such as formation of 
cluster, dipoles, quadrupoles, etc. This configurations 
can be reintroduced into the theory by allowing for the 
existence of equilibrium between monopoles and clusters. 
At the simplest level that we are considering in this pa- 
per, salt dipoles are taken as neutral species that do not 
interact electrostatically with PA, and can be ignored. 



APPENDIX B: THE TEMPERATURE T; 

In this Appendix we study the collapsed transition 
present at high temperatures and high concentration of 
salt. We also compare our results with those of Higgs and 
Joanny ||]. At high temperatures, the density of salt ps 
inside the volume Vi is slightly higher than po (see Fig.^) 
and, consequently, we can write 



Ps (1 + 5)po 



(Bl) 



where (5 <C 1. 

The density of free charges along the polymer can al- 
ways be defined as p = 7po. In this case, the total density 
of charges inside Vi is 



Pi = Po (1 + ^ + 7) 



(B2) 



where 7^/7. 

If Po < p, 7 > 1 and, therefore, at high tempera- 
tures, where association becomes irrelevant, the electro- 
static contribution to the free energy scales as p^^^. This 
attractive interaction is balanced by second virial term 
Wip^ /2, which forces PA to collapse. 

For high concentration of salt, 7 < 1. Therefore, in 
this case, 5 and 7 can be treated as small parameters. 
Solving Eq.(p3|) for 5 up to second order in 7 gives 



(5 = 7 



1 



AT* (zo + 1)^ 
1 



1 



1 



(1 - Bpo)2 4r* (zo -f 1)2 



16T* (zo + 1)4(1 - Spo)3 [(1 - Bpo)2 

n 3 r / 



1 



1 



zo 



AT* (zo + l)2 
1 



(3zo + l)(zo + 1) 



l-Bpo 



2T* 



zo(l - 3Bpo) 



(B3) 



Assuming that both the monomers and the ions of salt 
are point particles (B = 0) and expanding in zq, to lead- 
ing order the expression above reduces to Eq.(3.7) of 
Higgs and Joanny's paper Q. Unlike their analysis, how- 
ever, we find that when the full DH theory is used, the 
unphysical divergence at 1/(1 — i3po)^ = zo/4T*(zo -|- 1)^ 
is actually located inside the coexistence curve. Since the 
density of salt is always outside the coexistence region, 
this divergence is never encountered. 

The scaling form of th e radius of gyration is then given 
by substituting Eq.(p^) into Eq.(|22) 



3(a2 



1) 



W2P^ 



2(T*)2zo(zo + 1)2(1 -i3po)2 
1 1 Zo 



(1 - Bpo)2 AT* (zo + 1) 



p2 = . (B4) 



Note that this expression is just the minimization of the 
free energy of a usual polymer with the second virial co- 
efficient renormalized to 



Wi 



2{T*Y zo(zo + 1)2(1 -Bpo)2 
1 1 20 



(1 - Bpo)2 AT* (zo + 1)2 



(B5) 
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As in conventional polymers, for Wi > (good solvent 
case), the polymer will be fully extended, while for Wi < 
(bad solvent case) will assume a globular form. At 
T* = Tg, given by Wi = in equation above (theta 
solvent), there is a collapse transition, as shown in Fig.^. 



transition between the dilute globular state and the ex- 
tended state. Unfortunately, numerical analysis indicates 
that this transition happens at very low temperature. 



APPENDIX C: THE FIRST-ORDER PHASE 
TRANSITION AT LOW TEMPERATURES 
WHEN po < PC 

If the density of salt is chosen to be on the vapor side of 
the coexistence curve, at very low temperatures Eqs.(p2[), 
( p3| ) and (^^ exhibit two possible solutions: one with 
Ps = p" (low density of salt) and another with — 
(high density of salt) as T* 0. The first solution exists 
only for low temperatures, while the second is present 
for all temperatures. Consequently, we have to investi- 
gate if there would be a phase transition between these 
two phases as the temperature is decreased. In order to 
address this question, we have to compare the respective 
free energies for the two branches. 

The free energy associated with the solution with low 
density of salt is obtained by substituting 

x«l + 0(e-i/^^*), 

a«a + 0(e-i/2^*) , 
p,«p"+0(e-V2T*)^ 



(CI) 



into the Eq.(p^) what leads to 



Nf Nf 
'2T* 2^ttpmT* 



(C2) 



On the other hand, the free energy associated with the 
high density state is obtained by substituting 



a; « 1 + 



,-3/4^7rpMT*-lnT*/87rpj 



a « a + 0(1 - x) , 

Ps^P^ + 0{l~x) , (C3) 
into Eq. (|l9|) . The resulting free energy is given by 
.c ^ Nfx Nf 



2T* 



2\PkpmT* 



-O I ^ . (C4) 



Comparing Eqs.(|C^) and (| 

(3Q^ + (1 - x)fNj2T*, where 
exp{-3/4\/7rpMT 



InT^ 



gives 
(1 - 
/8t:pm}/VT^- 



that pg^ « 

x) scales as 
Since the last 



term in the above expression is positive for low temper- 
atures, the phase of low density of salt is more stable at 
low temperatures. This indicates that there is a phase 
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